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Some remarks on the calculation of the cohesive energy of 
ionic crystals 


By 5. O. Lunpevist and P. GO. FROMAN 


This paper is a complement to the earlier quantum mechanical treatments 
of ionic crystals in the Heitler-London approximation, which have been given 
by Hyutieraas, LANpsHorr and Lowprn. The method, which is restricted to 
crystals consisting of ions with complete electron shells, was founded by Hyt- 
LERAAS! in his investigation of LiH. He formulated the appropriate physical 
assumptions, the description of the crystal by a wave function having the cor- 
rect symmetry properties and the method of subtracting the energy of the free 
ions. LANDsSHOFF? has treated NaCl by using approximately ortho-normalized 
one-electron wave functions obtained as linear combinations of the one-electron 
orbitals for the free ions. He introduced the overlap integrals S,, defined by 
formula (2) below and ortho-normalized the wave functions correctly to the 
first order in the S,,, taking also certain second order terms into account. 
Later Lo6wp1n* has used a mathematically exact method for ortho-normalizing 
by means of the matrix calculus. However, in his calculations he used a series 
expansion in the S,, and considered an approximation similar to Landshoff’s, 
taking into account some further terms. Applying Lowdin’s formula for the 
cohesive energy without introducing this expansion, we shall show how the 
calculations can be performed even in those cases where the series expansion 
mentioned is not convergent. In this connection we have found it necessary to 
investigate the energy expression further. 

From the normalized wave functions wy, associated with the free ions, a set 
of ortho-normalized functions y, is obtained, according to Lowdin, by 


Qu = >ypal(1 +S) Haw. (1) 


The y, and g, are one-electron wave functions that depend on both space 
coordinates, r= (x,y,z), and spin coordinate, ¢. The transformation (1) is 
uniquely determined except for an arbitrary unitary transformation. S is the 
following Hermitian matrix introduced in this problem by Landshoff: 


1H. A. Hyuimraas, Z. f. Phys. 63, 771 (1930). 

* R. Lanpsnorr, Z. f. Phys. 102, 201 (1936); Phys. Rev. 52, 246 (1937). 

3 P. O. Lowpty, Dissertation, Uppsala 1948, Cf. also Arkiv f. mat., astr. o. fysik 35 A, 
N:o 9 and N:o 30; J. Chem. Phys. 18, 365 (1950). 
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Sit oF Syn) y (1) dt — On» (2) 


the sign fdt,:-- implying summation over the spin coordinate and integration. 
over the space coordinates.t Lowprn has supposed that the symbol (1 + S)7 i 
has a meaning. It should be remarked that for this to be true it is necessary | 
and sufficient that det (1 +S) == 0. We shall suppose this condition to be ful- | 
filled. In terms of the functions g, the energy is given by a well-known ex-_ 
pression. Substituting (1) into this expression, subtracting the energy of the | 
free ions and using the fact that the y, satisfy Hartree-Fock equations, Low- 

din finds the following formula for the cohesive energy”: 


where 
c a Jd, i el wi ( 
B= 5D D0 SH HOLM) + DD? er |G] wo) (3 a) 
g gy’ 99 bu le y 
BE, =— >) D>” (ur |@loy) (3 b) 
be v 


a 
| 


8 =— Dy Dy Purl |H| wy) +2 > vx |Elux)—(x |G] xu} (3) 
lb v zx 


B= DDD DPon Per (@Al@l ur) — AG] yah (3d) 


le v % 


The new symbols here have the following meaning: 


(| H|>) = fy (—2 3° —*— yds 


9 Ir: —ro| 


e 


(wv|G| 2A) =f f wi (1) wr (2) 5 yx (1) ya (2) dy dt, 


2|r;—r| 


where Z, is the charge of the nucleus g, and e is the electron charge. Pee 


. . . . . es 
denotes a summation over all indices vy except those associated with the ion 
with which w is associated. The meaning of oe as well as Ny (used later) 


hae d a ' 
is quite analogous. The quantities P,,, are the elements of a matrix P defined by 


S 
ie ) 


1 = fe 5 } * . . ‘ . 
; The symbol 1 in Yu (1) is written as an abbreviation for the position vector ry and the 
epi coordinate ¢, of an electron. 

We have written Ky, By and Hs, -- Ky instead of Léwdin’s Helstat, Hexch and Hs in order 
to indicate clearly that it is only the sum Ecoh=H, +B. +E3 + Ei, that has a real physical 
interpretation. 

* We find the notati ye) simpler i 

ation >) "" simpler than that introduced by Landshoff and taken over 


by Lowdin. 5 
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For the sake of simplicity we now specialize the treatment to LiH. This 
specialization is not essential, but is made only to get somewhat simplified for- 
mule. The symbol w has now reference to a definite lattice point g and a 
definite spin state j. We indicate this by writing m = (g,)), » = (g',9’) and so 
on. The given one-electron wave functions are written 


Pu(L) = pulri, $1) = YolTr) OF. (5) 
The quantities 


=f volr) po (r) dr — bgq 
form a matrix called s, the diagonal elements of which are evidently 0. We 
now define a new matrix p by 


Ss 
ies 


PS (6) 


which is possible if det (1 +s) +0. (Note that det (1 + S) = (det (1+ s))?.) 
The wy being real for LiH, the matrices s and p are both ayumerrcnt The 
matrices S and P are obtained by the formule 


Suv = Sqq' 03’ (7) 
Pus — Pagt 053"; (8) 
where  =(g,j) and v= (g’,j’). Using (3), (3a, b, c,d), (5) and (8), we find 


after some simplification the following expression for the cohesive energy of a 
Lil crystal!: 


Eeon = Ey + Eg + Ey + Ei, (9) 
where 
g 1 Z Z, , ¢ 13 , 
E,=2 > >” f ——— = 227 (99 9) + A\e9 go’) | (9 a) 
ging! ANG J 
2S ,) ( font 
E, = —@ > Dd (o9'9'9) (9b) 
Tend. 


Ss Pag’ a (9) {9 Za i(g' gg) pod 4 (g' ggg’ ‘ att 9 2 (g" g''g''g)} (9 c) 
g a 
\' \ ‘ LARUE ry eer ee, 3 
By= & >) 2) DD Paw" Poa 12(9"' 9" 99) — 9" 9'" Gg} (9 d) 
Pn 


g” og! 


The symbols appearing here are defined by 


1 The expression for an arbitrary ionic crystal consisting of ions with complete electron 
shells is obtained from (9), (9a, b,c, d) by replacing g by (g, 2) where n is an index enu- 
merating the electron orbitals of the ion g. 
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1 


ggg") = J ya (r) ie avy | yo (r)dr - 


Coke KE. 


ee ? L 
ggg’ 9") = SS pa (11) pa" (72) Lema} oat oh adele coms Ae 


From these definitions it is evident that 


(999) = (9 9'9)" 
and 


, VA ee HE Vf SIS: 


GV=G@ 0 Gey =6. 6.99) 


gg 99" )= (999 é: 

In order. to make full use of (9), (9a, b, c¢, d) we find it natural to consider 
the matrix p as a primary quantity, the calculation of which forms the firs’ 
step in the solution of the problem. 

If the condition? 


DAG | aa (10) 


g’ 


is fulfilled for all g, the matrix p can be expanded in the power series® thus, 


Peo Aes Cae (11) 


The corresponding expansion P = S — S* + §$?—+.--- mentioned in the intro- 
duction was introduced by Lownrn directly in (3c, d). We have performed some 
preliminary calculations on LiH and found that the series (11) is not conver- 
gent because of the large overlapping of the H-ions. Therefore we have tried 
to calculate p by another method, which seems to be generally useful in cal- 
culations according to the Heitler-London method. 

We write (6) in the following form: 


ee (Ogg + Sqq’) Dg'g' = Sg". (12) 


g 


This is a system of N? linear equations with N* unknown quantities pyy. (N is 
the number of ions in the crystal.) The matrices s and p are both sym- 
metrical, i.e., Sg'y = Sgy' and y'g = Pyg'. When g and g’ denote ions lying in 
the interior of a large crystal, the quantities py’ must also have symmetry 
properties corresponding to the lattice symmetry. This is evidently true if the 
perturbation from the surface layer of the crystal can be neglected, which ought 
to be the case if the crystal is sufficiently large. (Throughout this paper we 
are disregarding the influence of surface effects.) The coefficients dgq' + syg’ in 


* Evidently these relations remain true if the g are again interpreted as in the pre- 
ceding note. 
; As we consider a crystal of finite size, the matrix s is also finite. 
After we had arrived at the sufficient condition (10) for the validity of (11), Dr. Léw- 


din noticed and kindly pointed out to us that it can be found in Courant-Hilbert, Math. 
Physik I, p. 16. 
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(12) tend rapidly towards zero with increasing 7yy', and it seems probable that 
the same is also true for py’. After these arguments it seems promising to 
reduce the equation system (12) by making use of all symmetry properties of 
s and p. This reduced system can then be approximated by taking into ac- 
count only quantities associated with a few adjacent ions, whereby the number 
of resulting equations can be made rather small. We have tried this method 
and found it quite useful. The numerical calculations were carried out with the 
number of resulting equations ranging from 3 to 11. We found a satisfying 
convergence of the calculated matrix p with increasing number of equations. 
Further, the numerical results indicated that py, tends towards zero approxi- 
mately as rapidly as sy with increasing 1g’. 

As previous authors have done, we shall separate Hoon (9) into elementary 
parts, going, however, into more detail (especially as regards Hy) and making 
no neglections. F, (9a) can be written 


‘og Tog Tog 


ey Eg Eg! (gg! r aay ; 
By = 5D DOB 4 Mee) 4 7, Oeltee) _ a4 ai’yof> (9a) 
ee ae 
where e, = Z,— Ny is the valency of the ion g, and w,(r) is defined by? 
wy (7) = 82 f lyr (e) Pa? de —r f |yy(e) Pada}. 


This function is rapidly decreasing towards zero with increasing r. The new 
symbol in (9 a’) is defined by 


(00 0")o = fyb (re TED aay de 


For E£, we keep the formula (9b), whereas H; (9c) is transformed to the fol- 
lowing form: 


Eq’ Wg \?qq' , Poca’ 
E,= 2 DD” Dw West = &y! eo) + 2(99 go +2(999 0) + 
gg 99 


‘og 


¢ ~~ V(g) \V(g) Put ao ey wy , 
+22 >) D>" >)” paw fe" (9' 9" 9) * (99 Quo Fg Og G)}* (9c) 
g g gy" 


E, (9a’), E,(9b) and H; (9c’) are now written in such forms that every single 
sum in them is proportional to the crystal volume (if the lattice constant is 
kept fixed). 

The corresponding transformation of #,(9d) is more complicated, and there- 
fore a special investigation is motivated. We begin by considering the following 
parts of the sum 2e? 3 Ps py oe Dog" Do'g"'(9 g gg’): 


Tee! ae 


‘ For erystals other than LiH this definition of g(r) must be slightly generalized. 
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1°. The double sum D, that is obtained when g +9,9 =g9 and g”’ =g. 
To estimate it we note that for large distances rg’ we have 


(9900 Sis 


‘og 


Hence 
7 ge Tog 


D,=2¢ » Bait Pog Pg'g' (gg 99') wae 6 poral Poo Poe’ . 
gag 


5 . , ty TPS: g 
2°. The triple sum D, that is obtained when g +9,g° =g and g +{2 


To estimate it we use (6) and get (for large distances ry’) 


aes (99 ad) ~ DY? pyar S09" _ (P8)'o' — Po'o8oo' _ 
g Tgg' Tqg' 
_(s — P)a'a ~~ Po’ 9 899 Po'g' + Po'g Sag’ Poo’ 
Tog! Tog Tqg' 
Hence 
(g) (97) MD 9 2 IC) ye 
= 2 OO” doa Po'a (99" 99’) ~ —2e? >, Dd)” Pag Po'a’ 
Gino op ie tne Taq! 


ver a 


3°. The triple sum Ds, that is obtained when g’ + g, 9’ += ie and g° =@. 


By similar estimations as in 2° we get 
5) (99') , (9) ete 
Ds = 22 SO ner Oo 6g) = 28 DD) ate 
OG gy" Ta 5G Tog! 


Tog 


4°. The quadruple sum D, that is obtained when all four summation indices 
9,9, 9,9 are restricted to be different from one another. By using 
similar estimations as before, we get 


Di= = 2 e y Dae De ae Sor g' ) Dag" Da’ ¢ q” (g” ie ‘ag ")~ 2 e >) ae Pog Pa'g' : 


Sif g" qi oi) 9’ 


Tgg' 

The above estimations of D,, D., D; and Dy can be made more exact, but 
for the sake of brevity we omit the calculations required. For a discussion of 
the Di, we make the following remarks: It can be shown that, for fixed lattice 
constant, the sum 

pay (9) Pag Pg'g' (13) 
vz Tqg' 


re Vag 3 
* The symbol ae means a summation over all g” except g”’ =g and g”=g’. The 


” 


extension of this notation to the case when there are more restrictions is immediate. 
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tends to + co with increasing crystal volume V, unless the condition! 
Print + Poy = 9 (14) 
is satisfied. IH, on the other hand, (14) is satisfied, (13) is a sum of the Ma- 


delung type and has a finite limiting value. Our numerical calculation of the 
matrix p showed that (14) is not even approximately fulfilled (for LiH). Hence 


Di dDs ed. 
each one of the quantities ai a ae tends towards +°co or —oo with 
¢ ae 
increasing crystal size. However, their sum 7 Bs Yi ; cassia tends to 


a finite limit when V — oo. 
It is convenient to write the sum Hy = D, + D, + Dz; + Dy as follows: 


, : 4 f 1 
Ey=2e Dy 3 Pog Po' g' loo UIC G ow >} ne 
ig Mg 2 Too 
+4é 2 ae Pog jest Poo’ (99 99’) + me | - 
arian T9g' | 


+ 


y e > pe oar. aed Dag" Dy'g'” (g” gg’) ey ZeEG : (9 d’) 
gag g ge 


‘gy' . 
The rest of EH, (9d) can be written? 


Ey = E,—Ey=@ ee Pag (9999) + aa ie (499 Poo (9999) + 


g ip ee 
a9 (3(99 99) —(99'99')] — Pag Do'g' (9.9 9' 9)} 


( ( , (ae 
Di? DY? (Poa Poa" [3 (a9' 9" 9) — (a9 99") — 


oe Vine arp 

ee, .) @ NGI) NGS ye Wee A " 

— Pqg Pag (99' 9 9)} —@ >) >" path . > 0) yg" Dy gt (GG 99). (9.4) 
9 ig’ g’’ on 


The right-hand sides of (9d’) and (9d”) are now written so that every single 
sum in them is proportional to the size of the crystal. 

Previous authors? have classified the quantities according to their order of 
magnitude, essentially in the “overlap” of the atomic wave functions, and have 
neglected those of orders higher than an established one. This investigation 
shows, however, that only the complete expression for the cohesive energy has 
a real physical meaning and that great caution must be exercised in dividing 


} Tn the interior of a large LiH crystal the value of Pyg depends only on whether g is a 
Li’- or a H -ion. In the former case the value is denoted by Pyizi and in the latter 
by Pou: 

* The expression (9d) has been simplified by observing that the Li’- and H -wave func- 


tions Wg are real. 
® See p. 431, footnotes 1—3. 
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it into partial sums. Because of the agreement between the calculated and the — 
measured value of the cohesive energy for the crystals treated by Landshoff — 
and Léwdin it seems probable that the terms they neglected are relatively — 


small. For LiH the extension of the H-ions complicates the problem, and a 
new investigation of this crystal, formerly treated by Hylleraas, has been un- 
dertaken. The numerical results will be published later. 

Lastly we will make a remark about Hylleraas’ method of subtracting the 
energy of the free ions. If this energy is calculated by using the Hartree-Fock 


one-electron wave functions, which neglect ‘‘correlation”, we introduce an error | 
which can not be neglected when compared with the quantum mechanical 


part of the cohesive energy (i.e., the cohesive energy with the Madelung part. 
disregarded). We have verified this assertion only for LiH, but it ought to be 
general. Using the subtraction method of Hylleraas, we might expect that the 
energy error for the crystal and the energy error for the free ions cancel each 
other to a large degree. 
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